o 

(N 

o 

Q 



O 



^ 



X 



A New Microlensing Event in the Doubly-Imaged Quasar 

Q 0957+561 

Laura J. Hainline^, Christopher W. Morgan^, Joseph N. Beach^, C. S. Kochanek^, Hugh C 
Harris^, Trudy Tilleman^, Ross Fadely^, Emiho E. Falco^, and Truong X. Le^ 

ABSTRACT 

We present evidence for ultraviolet/optical microlensing in the gravitationally 
lensed quasar Q 0957+561. We combine new measurements from our optical 
monitoring campaign at the United States Naval Observatory, Flagstaff (USNO) 



r) • with measurements from the literature and find that the time-delay-corrected r 



band flux ratio rriA — rns has increased by ~ 0.1 magnitudes over a period of 



Ph' five years beginning in the fall of 2005. We apply our Monte Carlo microlensing 

O ■ analysis procedure to the composite light curves, obtaining a measurement of 

'^ • the optical accretion disk size, log{(rs/cm)[cos(i)/0.5]"'^/^} = 16.2 ± 0.5, that is 

, ^, ' consistent with the quasar accretion disk size - black hole mass relation. 

1> ' Subject headings: gravitational lensing: strong — gravitational lensing: micro 

QQ I accretion disks — quasars: individual (Q 0957+561) 

a^' 

p 

O ■ 1. INTRODUCTION 



The history of microlensing in the Zg = 1.405 quasar Q 0957+561 (hereafter Q0957), 
the first confirmed gravitational lens system, is paradoxically simple yet complex. The 
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possibilit y of microlensing in thi s double- image, wide-separ ation (~ 6^0 quasar was pro- 



posed by I Chang fc Refsdall ( 119791 ) soon after its discovery by IWalsh. CarsweU. &: Weymann 
( 19791). Controversy oyer the time delay, currently accepted to be Ajar = 417 days (e.g., 
Vanderriest et al.lll989l : ISchildlll990l : iKundic et al.lll997l : IShalyapin et al.ll2008l ). spawned nu- 
merous time variability studies of Q0957 at optical and radio wavelengths which produced 
a wealth of monitoring data. Microlensing of the quasar on long time scales was soon es- 
tablished by a change in the time delay-corrected magnitude diff erence of the two i mages 
of A{mA — fns) ^ 0.25 over t he co urse of five years, observed by I Vanderriest et al.l (119891 ) 
and confirmed by iPelt et al.l ( 119981 ). However, after that initial event, occurring between 
1982-1986, the difference light curve (the difference between the light curve for image A and 
the time-delay shifted light curve for image B) became nearly constant. Debate over the 
existence of measurable variability in the flux ratio percolated through the literature, with 
some authors reporting variability in the flux ratio on the order of several hundredths of 
a magnitude over timescales of weeks, days, and hours, while other groups simultaneously 
conc luded that no microlensing with amplitude \A(mA — wzr)! > 0.05 was detected 



see, 



e.g., iGil- Merino et al.l l2001t ISchmidt fc WambsganssI l2010l . for summaries and references). 
As the dispute surrounding the existence of short-timescale microlensing evolved, the ab- 
sence of long-time scale microlensing became increasingly evident from the difference light 
curve's steadiness over roughly 20 years. The only report of a possible long-time scale event 
after 1986, that of lOvaldsen et al.l (120031) of 5% - level y ariability over 300 days, was not cor- 
roborated independently b y IWambsganss et al.l (l2000l) with data taken over the same time 



(120081 ) find no evidence for microlensing in their 



period. Studies as recent as IShalyapin et al. 
data. 

With t he development of numerical simulation methods to analyze microlensing observa- 
tions (e.g., lKochanekll2004J : iBate. Webster, fc Wyithd 120071 : iBlackburne et al.l 1201 ll ). quasar 
microlensing has meanwhile become a powerful tool with which to quantitatively study the 
properties of lens galaxies and the structure of quasars. Although quasars and individual 
stars in lens galaxies cannot be resolved by conventional telescopes, by modeling the ampli- 
tudes of microlensing fluctuations we can measure the size of the continuum emission region, 
the masses and velo cities of the stellar microlenses, and the stellar-to-dark mass fraction in 



the lens galaxy (e.g.,lKochanekll2004j : IPooley et al.ll2007l : [Morgan et al.ll2008l : iMediavilla et al. 



20091 : iDai et al.ll2010l ). Through analysis of microlensing variability a t multiple wavelerigths. 



the surface brightness profile of the quasar may be measured (e.g., JPoindexter et al.l 12008 



Eigenbrod et al.ll2008l ). We are now moving into an era in which theoretical models of accre- 
tion disk structure can be tested with significant samples of lensed quasars, resulting in the 
observational confirmation of a relation between central black hole mass and accretion disk 



size ( JMorgan et al.ll2010l ). as well as the temperature profile of the lShakura &: SunyaevI (119731 ) 
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thin accretion disk model (e.g., iPoindexter et al.ll2008l : lAnguita et al.ll2008l : iMosquera et al. 
201lh . 



Since the central black hole in Q09 57 lies at the high end of the quasar black hole mass 
function (l.Ox 10^ Mg ^Assef et al.ll201ll ). Q0957 provides a potentially interesting and impor- 
tant data point for testing models of accretion disk structure. However, the available analyses 
of Q0957's only widely-acknowledged microlensing event from 1982-1986 pre-dated the use of 
the newer, highly quantitative microlensing analysis techniques, and the qualitative methods 
employed in the published analyses of the event produce either limits or constraints too broad 
to be informative for current tests of quasar s tructu re and studies of the surface mass den- 
sity of lens galaxies. For example, iPelt et al.l ( 119981 ) found that the 1982-1986 microlensing 
event was consistent with a quasar radius of "roughly" 3 x 10^^ cm and could be explained by 
microlens masses down to 10"^ Mq with only a small fraction (< 5%) in solar mass objects. 



Refsdal et al.l ( l2000l ) also analyzed the 1982-1986 microlensing event and subsequent 8-year 
quiescent period in Q0957 using three different simple lensing mass distribution models, ob- 
taining rough constraints on the range of microlens masses of 10^® < M/Mq < 5 and an 
upper limit on the size of the quasar's optical continuum region of R < 10^^ cm. Attempts 
to analyze the negligibly low microlensing variability in Q0957's difference light curv e with 
early numerical simulations (jSchmidt &: Wambsgansslll998l : IWambsganss et al.l 120001 ) could 
not constrain the quasar size. Such studies could only rule out microlens masses for a given 
quasar size, and in the end obtained results contradictory to analyses of the 1982-1986 
microlensing event. Thus, the lack of unambiguous microlensing events in Q0957 in the rel- 
atively recent era of computationally-intensive microlensing analysis has been unfortunate, 
preventing improved determinations of the quasar's size and microlens mass properties. 

However, as the saying goes, "good things come to those who wait." In this paper, 
we combine three new seasons of optical photometric monitoring of Q0957 with previously- 
published light curves to demonstrate the return of long-time scale uncorrelated variability 
in the light curves of Q0957. We take advantage of this new microlensing event to place 
the best constraints yet on the size of the optical accretion disk a nd the mea i i mas s of the 
microlenses using the Bayesian Monte Carlo analysis technique of iKochanekl (120041 ). In ^ 
and §2] we present our new monitoring observations of Q0957 and discuss the methods we 
use to model the uncorrelated variability in the light curves. We present our results in 
im Throughout our dis cussion we assume a flat cosmology with flu = 0.3, Qa = 0.7, and 
Ho = 70kms"^ Mpc"^ jHinshaw et al.ll2009h . 
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OBSERVATIONAL DATA 



We regularly monitor the flux of Q0957 A and B in the Sloan Digital Sky Survey r-band 
using the 1.55-ni Kaj Strand Astrometric Reflector at the United States Naval Observatory 
(USNO), Flagstaff Station, as part of the United States Naval Academy (USNA)/USNO 
Lensed Quasar Monitoring Program. Our program, which began in 2008, obtains three 
flve-minute exposures of the quasar per night, one night per week on average (weather 
permitting), using either the Tek2K CCD camera (O'/SSpixel"^) or the 2048 x 4096 EEV 
CCD camera (0'.'18 pixel" ). We present here data taken on 57 nights between 2008 March 
and 2011 June. The median stellar FWHM (seeing) of the images in our data set is I'.'S. 

We measure the quasar image fluxes using relative photome t ry. A detailed discussion 



of our image analysis methods was presented in iKochanek et al.l (|2006[ ). so we only briefly 
summarize our procedure here. We treat each quasar image as a point source and model 
the point spread function with three nested elliptical Gaussian components, keeping the 
relative astrometry fixed to that derived from Hubble Space Telescope (HST) images by 



Keeton et al.l ( l2000l ) (see Tabled]). The zi = 0.356 lens galaxy (JFalco et al.lll997l ) is modeled 
as a Gaussian approximation to a deVaucouleurs profile for which we fix the effective radius 
to that derived from the HST images of Q0957. The lens galaxy's flux is held constant in 
all epochs; we determine the optimal value for the constant flux by repeatedly fltting all the 
images as a function of lens galaxy flux and flnding the value which results in the lowest x^ 
sum over the entire data set. We compare the flux of each quasar image to four reference 
stars, located at (-59'.'2, -27"8), (-61'.'7, -109'.'1), (-68'.'5, -89'.'8), and (+111'.'8, -127^1) 
relative to Q0957 image A. Because these reference stars are signiflcantly redder than the 
quasar images {g — i ^ 0.8 — 0.9 for the stars and g — i ^ 0.2 for the quasar), and the 
two different detectors used in our USNO program have slightly different spectral responses 
across the r band, a small color term between the different detectors is introduced into our 
relative photometry measurements. We determined this color term by fltting the light curves 
of the quasar images obtained with the EEV camera to the light curves obtained from the 
Tek2K camera, yielding a magnitude offset of mTek2K — ^eev = 0.048 ± 0.002 mag. We 
apply this correction to the data from the EEV camera, and we include the contribution 
of the uncertainty in the color term, added in quadrature, in the photometric errors of 
measurements taken with the EEV camera. 

The measurements of Q0957 A and B from our USNO monitoring program are listed 
in Table [2l We note that we obtained data on two additional epochs which have not been 
included in Table |2l bringing the total number of observation epochs to 59. However, the 
images obtained on the two excluded dates were unusable due to high sky brightness levels 
on one night and corruption of the quasar images by cosmic rays on the other. 
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In order to obtain a longer time baseline for our microlensing analysis, we have supple- 
mented our new USNO light curves with published monitoring data. We use the i?-band pho- 
tometry from the 0.8-m telescope at the Instit uto de Astrofisica d e Cana r ias' (lAC) Te i de Ob - 
serv atory spanning the years 1996-2001 from ISerra-Ricart et al.l ( 119991 ) . lOscoz et al.l ( 120011 ) . 
and lOscoz et al.l (120021 ) , deriving the photom etric zeropoint offset by comparison to concur- 
rent r-band measurements of image A from iKundic et al.l (Il997l . Ams99-K97 = 0.04 mag); 
we neglect the small color term difference between r and R because in our microlensing 
analysis we use only t he relative flux ratio an d not the absolute flux. We also use the r-band 
monitoring data from IShalyapin et al.l ( 120081 ). spanning the years 2005-2007, and obtain the 
photometric offset (14.455 ± 0.018 mag) by analyzing their Liverpool Telescope/RATCam 
image^j for five different nights with our image analysis procedure and averaging the offsets 
in magnitude between the published photometry and the values produced by our pipeline. 
As the frequency of monitoring in those two data sets was nightly, on average, and our mon- 
itoring program obtains observations with a frequency closer to weekly, we have averaged 
the individual measurements into seven-day bins within individual observing seasons, after 
removing obvious (> 3cr) outliers in the data set, to prevent those data from having a dis- 
proportionate statistical weight relative to the new measurements from USNO. The center 
of each seven-day bin is defined as the mean HJD of the measurements included in that bin. 
To serve as errors on these seven-day averages, we compute the standard deviation of the 
individual data points included in each average. We show the final composite light curves 
for Q0957A and B from all of the data sources in Figure |T] Unfortunately, we are unable 
to find any published monitoring data for Q0957 which falls in the gap in Figure [T] between 
2001-2005; while the absence of data in this date range will affect our microlensing analysis, 
we will still derive meaningful results. 



We investigated including the two seasons of r-band monitoring from iKundic et al. 



( 119971 ) from 1995-1996 in our composite light curve for Q0957, which would have extended 
our time axis even further. However, we chose to not include this data in our analysis 
because the flux contribution to image B from the lens galaxy has not been subtracted from 
the photometry in that data set. Moreover, archival images from this data set are no longer 
available at the Apache Point Observatory or any of its operating consortium institutions, so 
it was not possible to carry out the deconvolution of image B from the lens galaxy ourselves. 



^ Shalyapin et al.l have made their monitoring images pubhcly available at the website 



http://dc.zah. uni-heidelberg.de/liverpool/res/rawframes/q/form 
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3. ANALYSIS 

3.1. Difference Light Curves and Microlensing Signal 

To analyze the light curves of Q0957 for the presence of uncorrelated variability, we 
must first eliminate the variability intrinsic to the quasar source itself. We accomplish this 
by first shifting the light curve of image B by the system's 417-day time delay and then 
performing a linear interpolation of image B's shifted light curve to generate a set of pho- 
tometric measurements at the same epochs of observation as those in image A's (unshifted) 
light curve. We discard any data points that were interpolated in the inter-season gaps. 
Finally, we subtract from the light curve of image A the shifted light curve of image B, 
creating a time-delay-shifted difference light curve in which only the uncorrelated variability 
remains o 

In Figure [2] we present the time-delay-shifted difference light curve for Q0957, focusing 
on the last six years of data. We note in advance that the frequency of observations of Q0957 
in the data sets we use does not p ermit us to confirm o r refute the reported observations of 
intra-day microlensing variability ( ICoUey fc Schildll2003l ). Moreover, the size of the error bars 
on the measurements render insignificant hints of fiux ratio variability on time scales of weeks 
or mon ths, so we cannot examine reports of the short time scal e microlensing variability in 
Q0957 JSchildlll996l : ICoUev fc Schildlboooi lOvaldsen et"aDl2003h . We do observe a slow but 
steady brightening of i mage B relative to image A in the r band, beginning at the start of the 
Shalyapin et al.l (120081 ) data near MH JD 3650 (in the frame of image A) and continuing to the 
limit of our USNO data set for image B at MHJD 5311. This clear signature of microlensing 
is evolving quite slowly, as indicated by the slope of a least-squares straight line fit to the 
difference light curve of 0.016 ± 0.006 magyr"^ (quoted uncertainty is 3 a). Such a slo w 



drift in fiux ratio would not have been statistically significant at the time IShalyapin et al. 



published their data, if even noticeable, as they could only constrain the fiux ratio up to date 
MHJD 3831 in the frame of image A with their data set; our monitoring data from USNO 
provide the extension in time necessary to confirm that the change in fiux ratio is significant 
in comparison to the observational uncertainties. The duration of the complete microlensing 
event is presently unclear from the shape of the difference light curve, but the nearly five 
year span (after time delay correction) so far implies the time scale will be at least that long. 
We are thus fortunate to be monitoring the first significant multi-year microlensing event in 
Q0957 in twenty years, and the first to be observed since the system's time delay was firmly 
established. 



^We note that if we instead shift and interpolate the hght curve of image A to the observational epochs 
of image B, no significant changes in the difference light curve or analysis results are observed. 
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3.2. Quantitative Monte Carlo Microlensing Analysis 

The long time scale and low amplitude of the new microlensing varia bility are qualita- 
tive ly consistent with th e relatively large optical quasar size predicted by iPelt et al.l (Il998[ ) 
and iRefsdal et al.l (120001 ) in their analyses of the 1982-1986 microlensing event and are not 
surprising given that the black hole is very massive. We can place new, improved quantitative 
constraints on the quasar size and the mass and velocity of the microlenses by analyzing the 
mi crolensing variab ility in our compiled light curve of Q0957 with the techniques described 



m iKochanekl (J2004J ). The analysis method has three major components. First, with the 
selection of a model for the macroscopic (strong) lensing, a stellar mass function to describe 
the microlens mass distribution, and an accretion disk model, we generate microlensing mag- 
nification patterns for a range of mass contributions from the dark matter halo of the lens 
galaxy. Next, we use a Monte Carlo method to generate large numbers of trial light curves 
from the magnification patterns for random source (quasar) trajectories, and fit the simu- 
lated light curves to the full observed light curves (years 1995-2011). Finally, we perform 
a Bayesian statistical analysis on the goodness-of-fit (x^) statistics of the light curve fits to 
calculate probability distributions for accretion disk size, effective source velocity, and mean 
microlens mass. 

To descr i be th e macroscopic lensing in Q0957, we utilize the results from the study of 
Fadely et al.l ( l2010l ). who develop detailed models of the quasar lens using as constraints faint 
knots and structures visible in quasar- and lens-subtracted HST images of the strong lensing 
region as well as the results of w eak lensing analysis. The lens potential in Q0957 has been no- 
toriously difficult to mode l (e.g.. lKochaneklll99ll : lGrogin fc Narayanlll996l : lBernstein fc Fischer 
19991 : iKeeton et al.l |2000|) , at least in part due to the two-ir aage nature o f the system and 
the presence of a galaxy cluster surrounding the lens galaxy. iFadely et al.l construct a series 



of models containing an isophotal model of the lens galaxy as the stellar component, a con- 
centric elliptical dark matter halo component representing the lens galaxy and the cluster 
halo, and a set of general third-order terms in the Taylor series expansion of the potential 
from the lens env ironment. Here, we use the convergence (k) and shear (7) of the model 
from lFadely et al.l which lies at the peak of the posterior probability distribution to generate 
the microlensing magnification patterns. So that we may marginalize over the uncertainty in 
the dark matter fraction in our analysis, we create a series of ten models in which the stellar 
component of the local convergence k^, in the vicinity of image B varies linearly in the range 
0.1 < {h^:/k)b < 1-0. Since image A is much farther from the lens galaxy, we constrain its 
stellar mass fraction (k^,/k)^ such that the ratio (fi;*/fi;)^/(/€*//€)_B is held fixe d to its value 
in the iFadely et al.l model (~ 0.2). We note that the flux ratio favored by the iFadely et al. 
models (fh/fa ^ 0.53) is somewhat lower than what is actually observed in the optical (1.1; 



Shalyapin et al.ll2008t ICoicoechea et al.ll2005l ). 



The magnification patterns are 4096 x 4096 images with an outer scale of 20{Re), where 
(Re) is the Einstein radius for the mean microlens mass (M) projected into the source plane. 
The outer dimension and pixel scale are chosen to be sufficiently large to representatively 
sample the magnification patterns and sufficiently small to resolve the accretion disk. The 
stellar mass function we use to generate the population of microlenses for the patterns is a 
power law, dN/dM oc M~^'^, with a ratio of maximum-to-minimum mass (dyn amic range) o f 
50. This function reasonably approximates the Galact ic disk niass fu nction of lGouldl (120001 ): 
assuming a different mass function (such as that of ISalpeten Il955l ) will have a negligible 



effect on our results due to the other, larger sources of uncertainty (see, e.g.. |Paczyhskilll986 



Wyithe et al.ll2000l ). We create four independent magnification patterns for each quasar 
image for all ten evenly spaced stellar mass fractions. 

Because the effects of finite continuum source size are generally not negligible for quasar 
microlensing, we convolve the magnification patterns with the surface brightness profile of 
the source for a grid of source sizes before computing trial light curves. We model the 
accretion disk as a face-on, thin disk radiating as a blackbody with a power-law temperature 
profile T oc i?^/"^. The scale radius of the disk, r^, is defined as the radius at which the 
disk temperature matches the rest-frame wavelength of the filter used in our monitoring 
observations, kT = he/ X^est (for r-band monitoring of Q0957, Aresi^ = 2593 A). Our model 



matches the outer regions of the thin disk model of IShakura &: SunyaevI (Il973l ). but we 
neglect the drop in temperature in the center due to the inner edge of the disk and the 
correction factor from general relativity to avoid introducing additional parameters. The 
effect of the simplification on our result for the accretion disk size is small compared to the 
uncertainties in other p arameters, as lo ng as the disk is significantly larger than the radius 
of the inner disk edge (JDai et al.l 120101 ) . We caution that because microlensing amplitude 
depends on the projected area of the source and not its shape, and because we assume a 
face-on disk model, the disk radius we infer will be an effective radius defined as the radius 
of a circle of the same projected area as the accretion disk. The true radius will be (cos i)~^/^ 
times this effective radius, where i is the inclination of the disk. 

We generate 10^ trial light curves for each of the 40 sets of magnification patterns. In 
each realization we randomly select an initial position and effective velocity for the source 
trajec tory, assuming the values are independent and uniformly distributed. As in iKochanek 
(120041 ). we neglect the motion of the stars within the lens galaxy and describe the observer's 
motion as the projection of the CMB dipole velocity onto the lens plane. We then use 
Bayesian methods to analyze the x^ statistics of the many light curve fits and calculate 
likelihood functions for the quasar source size and velocity in Einstein units (f^ and Ve) since 
all results scale with the unknown mass of an average microlens (M). We estimate the true, 
unsealed physical source size r^ = fs{M/MQY^'^ and mean microlens mass (M) by combining 
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the probability distributions of the scaled source size, P{rs), and the scaled effective velocity, 
P{ve), obtained from the Bayesian analysis with a prior probability function for the true, 
un sealed effe c tive s ource velocity, P{ve). We construct P(fe) using the method described 



m 



Kochanekl (120041 ). a pplying the measured stellar velocity dispersion for the lens galaxy 



(a^: = 288 ± 9kms ^; iTonry fc Franxl 119991 ). and obtaining the dispersion of the peculiar 
vel ocity distribution at th e redsh ifts of Q0957 and the lens galaxy from the power-law fits 



by iMosquera fc Kochanekl ( 1201 ll ) to the peculiar velocity models of Tinker et al. (2011, in 



preparation) . 



4. RESULTS AND DISCUSSION 

In Figure |3] we show the probability distribution for the scaled effective source plane 
(Einstein) velocity Ve resulting from the Bayesian analysis of our Monte Carlo light curve 
simulations, along with the probability density for the true source velocity, -P(fe), used to 
determine the mean microlens mass (M) and physical source size r^. The v^ distribution 
has a median of 1600 km s~ and is wide, with 68% confidence range of 600 km s^ < Ve < 
3500 km s~ . The distribution is also visibly asymmetric, with a substantial low- velocity tail. 
The width and asymmetry indicate the wide range of feasible solutions to the microlensing 
variability in Q0957 and reflect the ease with which the system's light curve was fit by the 
Monte Carlo code. Because the Monte Carlo simulations produce large numbers of acceptable 
fits and the x^ P^r degree of freedom for the best fits tends to be quite low (~ 0.4) even 
with insignificant (small) assumed values for the systematic errors of the flux ratio and 
photometry, we conclude that the spread in the velocity distribution is a consequence of the 
low intrinsic variability of Q0957 coupled with the very low amplitude of the microlensing 
signal, making the light curve simple to reproduce. An additional factor may be that the 
lack of observational constraints in the four-year gap in the compiled light curve permits 
considerable freedom in model behavior, as shown in Figure HI We tested this hypothesis 
by conducting an experiment in which we forced the flux ratio in the gap to stay within 
ibO.05 mag of its mean value in the dates bracketing the gap, (rra^ — m^) = 0.04 ± 0.03 mag. 
The experiment resulted in a minor narrowing of the Ve distribution, but the change was 
not significant enough to justify any firm conclusions about the influence of the gap on our 
results. 

Because the mean microlens mass is proportional to the inverse square of the scaled 
effective velocity {ve = fe(M/M0)~^/^), the broad range in permitted velocities causes the 
mean mass to be poorly constrained, which is apparent in the probability distribution for 
(M/Mq) shown in the right panel of Figure |3l The median of the (M/Mq) distribution is 
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0.2 Mq, with a 68% confidence range of 0.02 M© < (M) < 1.3 Mq. Our estimate of tlie mean 
microlens mass is appropriate for stars in an old e lliptical galaxy, and falls well within the 
most probable range found by iRefsdal et al.l (120001 ) of 10~^ < M/Mq < 5. Our constraints 
do not permit a dominant population of planetary and sub-plane tary mass (< 10~^ Mc?)), 
compa ct m icrolenses in the lens galaxy, consistent with the results of lSchmidt fc Wambsganss 
(Il998l ) a nd IWambsganss et al.l ( 120001 ). Our new constraints are also not consistent with the 
result of iPelt et al.l (119981 ) for the 1982-1986 event, in which a significant population of 
stellar mass microlenses is ruled out. However, as our microlensing analysis technique is 
significantly more sophisticated than that used in earlier studies and does not rely on the 
assumption of a single value for the source size or effective velocity to derive the microlens 
mass, we consider our new result to be much more robust. 

While we have attempted to constrain the stellar-to-dark mass fraction k^^/k in the 
lens galaxy of the Q0957 system, the probability distribution resulting from our Bayesian 
analysis is uninformative, with no strong peaks or trends favoring any particular k^^/k, value. 
We suspect that the primary reason for this failure is that there is insufficient uncorrelated 
variability in the light curves of Q0957 to constrain the dark matter fraction. However, the 
addition of X-ray data to the m icrolensing analysis may provide stron ger constraints, as it 
has in the case o f PG 1115-1-080 jMorgan et ahlbood : [Poolev et al.ll2009h and RXJ 1131-1231 
JPai et al.ll201o[ ). 



In Figure we show the probability distribution for the physical source size r^ of the 
quasar's accretion disk in the observed-frame r-band which results from our microlensing 
analysis. The distribution as shown has not been corrected for inclination i; however, in 
the text that follows all the numerical quantities we discuss will be corrected by a factor 
of (cosi)"^/^ assuming an inclination of 60°, the expectation value of a random distribution 
of inclinations. Note that the thin disk scale radius can be converted to a half-light radius 
using the relation ri/2 = 2.44 r^. The median of the microlensing thin-disk size distribution 
is log(rs/cm) = 16.2 ± 0.5, where the error bar represents the bounds of the 68% confidence 
interval. When converted to a h alf-light radius, which is comparable between different source 
models ( iMortonson et al.ll2005[ ). the source size we obtain from our analysis [log(ri/2/cm) = 
16.5 ± 0.5] appears to be marginally consiste nt with the i^-ban d half-light ra dii (converted 
from a of a Gaussian disk profile) obtained by lPelt et al.l ( Il998l . 10^^'^ cm) and lRefsdal et al. 



( I2OOOI . < 10^^'^ cm) from the data set covering the 1982-1986 microlensing event and the 
following 8 years of stabl e flux ratio. We note th a t our sour ce size result is substantia lly 
larger than those used by ISchmidt fc Wambsgansd (119981 ) and IWambsganss et al.l (I2OOOI ) to 
constrain the microlens mass distribution [log(rs/cm) < 14.2]; yet our mean microlens mass 
is still consistent with their constraints. We also experimented with the use of a uniform 
prior on the microlens mass of 0.1 < {M /Mq) < 1.0. Although the median of the resulting r^ 
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distribution is essentially the same as that found without imposing the mass prior, within the 
uncertainties, the distribution yielded by the prior contains a small but significant probability 
in solutions with small source sizes which are not physically reasonable. We suspect that 
the occurrence of this "shelf at the low end of the size distribution is an artifact of a 
combination of the paucity of microlensing variability in the observations and the mass 
prior itself. By limiting (M) to a range about which its distribution (derived from the Ve 
prior) is not symmetric, we give more statistical weight to solutions from one side of the 
distribution. Since the Einstein velocity (-Oe) scale implied by a given microlens mass scales 
as Ve = Ve{M/MQ)~^'^, the solutions with very low Einstein velocities are given more weight 
relative to the highest velocity solutions when the prior is applied. In the case of Q0957, the 
solutions with low Einstein velocity correspond preferentially to small source sizes; however, 
at very small sizes all solutions become equally likely since none of them pass near a caustic. 



Morgan et al.l ( l2010l ) used the same microlensing analysis technique as we use here on 



light curves of a sample of 11 different gravitationally lensed quasars to show that the quasar 
accretion disk size at 2500 A inferred from microlensing (r25oo) is correlated with central black 
hole mass as r25oo oc M^^'^'^'^'^ . Although the scaling of the observed correlation is consistent 
within the uncertainties with thin disk theory (R oc Mg(j), the correlation implies a quasar 
radiative efficiency rj that is approximately an order of magnitude lower than is expected 
based on the local supermassive black hole mass and quasar luminosity functions. We can use 
this r25oo^AfBH correlation as an independent check on our assertion that the uncorrelated 
variability we have observed in Q0957 is a result of microlensing. To scale the observed- 
frame r-band source size we find for Q0957 to a rest-frame wavelength of 2500 A, we assume 
the R oc A"^/^ scaling of thin disk theory which is supported by several observational studies 



of quasar microlensin g (e.g., lAnguita et al.l 120081 : iPoindexter et al.l 120081 : iFloyd et al.ll2009 
Mosquera et al.ll201ll ). In doing so, we note that the effect of the wavelength scaling is 



minimal for Q0957 since the effective wavelength of the observed-frame r filter corresponds 
to 2593 A at z = 1.41. For the black hole mass we use the CiV emission-line estimate of 



Assef et al.l (l201l[ ). which has an associated systematic uncer tainty of 0.33 d ex. In Figure |6l 
we place Q0957 on the best-fit r25oo^^BH relation found by iMorgan et al.l . along with the 
quasar data points used to derive the correlation. Q0957 is consistent with the relation; 
the agreement supports our assertion that microlensing is the source of the quasar's long- 
timescale uncorrelated variability. 



Similarly to iMorgan et al.l ( l2010l ) , we also show for comparison in Figure [6] quasar source 
sizes calculated by two additional methods. First, we calculate the theoretical scale radius of 
a thin accretion disk at a rest wave length of 2500 A from its central black hole mass ( "theory 
size") using Eq. 2 of IMorgan et al.l a r id ass uming accretion efficiency r; = 0.1 and Eddington 
ratio L/Le = 1/3 (JKoUmeier et al.ll2006l ). Second, we calculate the disk size under the 
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thin disk model (a blackbody with a T oc R~^^^ temperature profile) constrained by the 
observed magnification-co rrected optical flux/luminosity from HST F814W measurements 
("flux size", see Eq. 3 of [Morgan et al.l ). When we compare the results of the different 
methods of size calculation for Q0957, we observe that the micro l ensin g size is larger than the 
theory si ze, also noted for oth er lensed quasars by [Morgan et al.l (120101 ) and lBlackburne et al. 
(l201ll ). iPooley et al.l (120071 ) arrive at a similar conclusion, although they use bolometric 
luminosity-based black hole masses to calculate their theory sizes which are thus more similar 
to our flux sizes. Furthermore , the microlen sing size for Q0957 is l arger than its flux size as 
well, which was also found by iMorgan et al.l and iMediavilla et al.l ( 120111 ) for their respective 
quasar samples. However, the signiflcance of the disc repancy we flnd for Q0957 is not as 
large as that typically found by lBlackburne et al.l (120111 ) for their sample: while they rule out 
the black hole mass-based theoretical prediction by at least 3a in nearly every case, the flux 
size for Q0957 falls within the 68% confldence interval (essentially la) of the microlensing 
size, and the theory size within the 80% confldence interval (~ 1.3a). While this may simply 
reflect the relatively large uncertainty in the microlensing size estimate for Q0957 compared 
to other microlensed quasars, it is also possible that the difference in signiflcance is related 
to our different treatment of systematic errors from the analysis of Blackburne et al. The 
origin of th e discrepancies bet ween the microlen sing size, theory size, and flux size remain 
unclear. As lMorgan et al.l and lBlackburne et al.l point out, adjusting the mass accretion rate 
(L/tiLe) of Q0957 by lowering rj or using a higher-than-typical Eddington ratio may resolve 
the discrepancy of the theory size with the microlensing size, but such adjustments will not 
address the fact that the flux size is smaller than the microlensing size. 

One possible explanation for the size discrepancy is that the observed r-band flux may 
be contaminated by (a) UV/optical light from the continuum source scattered by the broad- 
line region (BLR); or (b) higher energy continuum emission reprocessed by the BLR and 
re-emitted as UV/optical emission lines. Contamination from line emission is an especially 
strong concern for our observations of Q0957 because the quasar's redshifted Mgll line falls 
within the r bandpass and could cause us to overestimate the microlensing size. Thus, it is 
important that we consider contamination scenarios in our analysis. So, we have repeated 
our Monte Carlo light curve simulations assuming that different percentages of the observed 
r-band flux can be attributed to unmicrolensed contamination from the BLR. When we do 
so, we flnd that if we assume 10% contamination we still obtain a median physical accretion 
disk size in r-band of log(rs^io/cm) = 16.1 ± 0.5. Even when we assume a 30% contribution 
to the observed flux from contamination, we obtain log(rs^3o/cm) = 16.0^q4, insigniflcantly 
different from our original microlensing size. Thus, neither scattering nor line contamination 
appear to produce a signiflcant effect on the microlensing size result, relative to our other 
sources of uncertainty. Moreover, even if line contamination produced a signiflcant decrease 
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in the micro lensing size, the flux size would be reduced along with the microlensing size, and 
the two would still not be reconciled. 

Although we are not able to resolve the discrepancy between the different accretion disk 
size estimates and our microlensing size measurement, we are nonetheless very encouraged 
to observe a long-term microlensing event in a quasar at the high end of the black hole mass 
function. Moreover, Q0957 is relatively easy to monitor with ground-based observatories, 
since the quasar images and lens galaxy are relatively widely separated and the time delay 
is well-determined. Future mult i- wavelength optical monitoring as well as X-ray monitoring 
will be most informative to examine the temperature profile of the accretion disk and compare 
the X-ray and optical sizes. 

The authors wish to thank L. Goicoechea, E. Turner, & R. McMillan for their assistance 
in obtaining historical photometric data for Q0957. We also thank the anonymous referee 
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C.W.M.), AST-0708082, and AST-1009756 (to C.S.K.). C.W.M. also gratefully acknowl- 
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Table 1. HST Astrometry and Photometry of Q 0957+561 



Component 



Astrometry 
AR.A. ADec. 

farcsec) (arcsec) 



F555W 

(mag) 



Photometry 

F814W 

(mag) 



F160W 

(mag) 











A 

B +1.229 ±0.005 -6.048 + 0.004 

Lens Galaxy +1.406 + 0.006 -5.027 + 0.005 



17.09 + 0.07 
17.11 + 0.05 
19.05 + 0.04 



16.71 + 0.05 
16.78 + 0.04 
17.12 + 0.02 



15.60 + 0.03 
15.68 + 0.03 
15.14 + 0.09 



Note. — HST astrometry and photometry for Q0957 are taken from lKeeton et al.l (J2000l ) 
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Table 2. Q 0957+561 Light Curves 



HJD - 2450000 


Seeing 


QSO A 


QSOB 


(Stars) 


(days) 


(arcsec) 


(mag) 


(mag) 


(mag) 


4554.712 


1.3 


2.671 ± 0.004 


2.419 ± 0.004 


0.013 ± 0.002 


4555.767 


1.4 


2.659 ± 0.004 


2.412 ± 0.004 


0.010 ± 0.002 


4561.722 


1.3 


2.665 ± 0.004 


2.417 ±0.004 


0.013 ±0.002 


4570.669 


1.5 


2.669 ± 0.004 


2.413 ± 0.004 


-0.004 ±0.002 


4584.677 


1.4 


2.664 ± 0.004 


2.378 ± 0.004 


0.008 ± 0.002 


4596.704 


1.4 


2.682 ± 0.004 


2.377 ± 0.005 


-0.002 ± 0.002 


4613.691 


1.4 


2.689 ± 0.004 


2.392 ± 0.004 


0.008 ± 0.002 


4617.702 


1.3 


2.691 ± 0.004 


2.412 ± 0.004 


0.010 ± 0.002 


4793.005 


0.9 


2.733 ± 0.004 


2.541 ± 0.005 


0.007 ±0.002 


4802.980 


1.4 


2.711 ±0.005 


2.544 ± 0.005 


-0.008 ± 0.002 


4807.929 


1.3 


2.728 ± 0.005 


2.563 ± 0.005 


-0.002 ± 0.002 


4829.857 


1.3 


2.747 ± 0.004 


2.604 ± 0.004 


0.007 ±0.002 


4833.894 


1.4 


2.726 ± 0.004 


2.603 ± 0.004 


0.006 ± 0.002 


4839.966 


1.6 


2.722 ± 0.005 


2.608 ± 0.005 


-0.017 ±0.002 


4846.758 


1.8 


2.730 ± 0.004 


2.621 ± 0.004 


-0.011 ±0.002 


4862.882 


1.2 


2.706 ± 0.004 


2.631 ± 0.004 


0.013 ±0.002 


4877.716 


2.0 


2.684 ± 0.004 


2.618 ± 0.004 


-0.009 ± 0.002 


4883.932 


1.1 


2.673 ± 0.002 


2.624 ± 0.002 


-0.078 ± 0.002 


4891.824 


1.3 


2.635 ± 0.004 


2.625 ± 0.004 


0.054 ±0.001 


4911.800 


1.1 


2.618 ±0.004 


2.665 ± 0.004 


0.013 ± 0.002 


4942.797 


1.7 


2.573 ± 0.003 


2.655 ± 0.004 


0.011 ±0.001 


4964.715 


1.2 


2.587 ± 0.004 


2.630 ± 0.004 


0.013 ±0.001 


4981.695 


2.2 


2.580 ± 0.005 


2.608 ± 0.005 


-0.028 ±0.002 


4997.677 


1.2 


2.608 ± 0.004 


2.626 ± 0.004 


0.013 ±0.002 


5144.981 


1.3 


2.573 ± 0.004 


2.656 ± 0.004 


0.004 ± 0.002 


5157.964 


1.1 


2.578 ± 0.004 


2.642 ± 0.004 


0.001 ± 0.002 


5186.024 


1.2 


2.558 ± 0.004 


2.603 ± 0.004 


0.005 ± 0.002 


5201.946 


1.2 


2.579 ± 0.004 


2.639 ± 0.004 


0.003 ± 0.002 


5208.993 


1.2 


2.585 ± 0.004 


2.652 ± 0.004 


0.009 ± 0.002 


5241.874 


2.1 


2.627 ± 0.004 


2.668 ± 0.004 


-0.010 ±0.002 


5261.687 


1.4 


2.631 ± 0.004 


2.667 ±0.005 


-0.016 ±0.002 


5296.691 


0.8 


2.641 ± 0.004 


2.617 ±0.004 


0.015 ±0.002 
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Table 2 — Continued 



HJD - 2450000 Seeing QSO A QSO B (Stars) 

(days) (arcsec) (mag) (mag) (mag) 



5311.690 


1.3 


2.634 ± 0.004 


2.575 ± 0.005 


-0.011 ±0.002 


5320.745 


0.9 


2.639 ± 0.004 


2.546 ± 0.004 


0.017 ±0.001 


5324.729 


0.7 


2.636 ± 0.004 


2.536 ± 0.004 


0.018 ±0.001 


5332.717 


1.1 


2.630 ± 0.004 


2.528 ± 0.004 


0.065 ± 0.002 


5348.695 


1.1 


2.605 ± 0.004 


2.480 ± 0.004 


0.013 ±0.001 


5352.685 


0.9 


2.596 ± 0.004 


2.483 ± 0.004 


0.014 ±0.002 


5358.688 


1.8 


2.560 ± 0.004 


2.463 ± 0.005 


-0.003 ± 0.002 


5506.991 


1.4 


2.650 ± 0.004 


2.474 ± 0.004 


0.007 ±0.002 


5519.030 


0.9 


2.712 ± 0.003 


2.490 ± 0.003 


-0.008 ± 0.001 


5527.937 


1.7 


2.692 ± 0.004 


2.483 ± 0.004 


-0.008 ± 0.002 


5535.032 


1.1 


2.704 ± 0.006 


2.489 ± 0.006 


-0.016 ±0.002 


5543.989 


1.0 


2.713 ± 0.004 


2.497 ± 0.004 


0.014 ±0.002 


5563.958 


1.3 


2.705 ± 0.004 


2.477 ± 0.004 


0.009 ± 0.001 


5590.779 


1.0 


2.729 ± 0.004 


2.496 ± 0.004 


0.017 ±0.001 


5604.847 


1.5 


2.711 ±0.004 


2.465 ± 0.004 


0.004 ± 0.002 


5621.687 


1.0 


2.705 ± 0.004 


2.482 ± 0.004 


0.012 ± 0.002 


5626.779 


1.2 


2.698 ± 0.005 


2.487 ±0.005 


0.003 ± 0.002 


5649.715 


1.1 


2.661 ± 0.004 


2.538 ± 0.004 


0.015 ±0.001 


5653.710 


1.1 


2.659 ± 0.004 


2.545 ± 0.004 


0.014 ±0.001 


5664.750 


1.2 


2.658 ± 0.004 


2.559 ± 0.004 


0.003 ± 0.002 


5674.748 


1.6 


2.648 ± 0.008 


2.566 ± 0.008 


-0.091 ± 0.003 


5684.672 


1.3 


2.653 ± 0.004 


2.564 ± 0.004 


-0.007 ±0.001 


5702.677 


1.7 


2.663 ± 0.004 


2.572 ± 0.004 


0.005 ± 0.002 


5712.673 


1.0 


2.671 ± 0.004 


2.582 ± 0.004 


0.012 ±0.002 


5733.651 


1.2 


2.673 ± 0.005 


2.566 ± 0.005 


-0.065 ± 0.002 



Note. — HJD is the Heliocentric Julian Day. The magnitudes listed in the 
QSO A and B columns are measured relative to the comparison stars. The 
magnitudes in the (Stars) column are the mean magnitudes of the comparison 
stars for that epoch relative to their mean over all epochs. 
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Fig. 1. — Compos i te r-b and light cu rve of Q0957 including; h istorical me asurements from 
Serra-Ricart et all Jl999l . lAC data), lOscoz et all J200ll . bood lAC data), Ishalvapin eTaL 
(I2OO8I . Liverpool data), and our new USNO data. The ~ 1600 day gap in the center of the 



light curve reflects an absence of published photometric monitoring for this system, but will 
not prevent us from deriving useful constraints from our microlensing analysis. 
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Fig. 2. — Time-shifted difference light curve of Q0957 from 2005-2010, displaying the slow 
increase in flux ratio of images A and B over five years that is a clear signature of microlensing. 
A simple linear fit to the observed data points indicates an average change of ~ 0.02 magyr"^. 
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Fig. 3. — Left panel: Probability distribution for the effective source velocity Ve for Q0957. 
The heavy line is the scaled effective velocity distribution in Einstein units, which has median 
Ve = 1600 km s~ , while the light line indicates the prior probability distribution for the true 
source velocity. Right panel: Probability distribution for the mean microlens mass (M) 
for the lens galaxy in the Q0957 system, which has a median value of (M) = 0.2 M0. 
(M) is calculated by combining the scaled effective source velocity distribution with the 
prior probability distribution for the true source velocity. The width of the (M) distribution 
reflects the uncertainty in the effective velocity apparent in the left panel, since (M) depends 
on the inverse square of the scaled effective velocity -Og. 
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Fig. 4. — Time-shifted difference light curve of Q0957 with a variety of simulated light curves 
that are good fits to the observed data, illustrating the ease with which the observed data 
are fit with simulated microlensing trajectories for a variety of different physical parameters. 
The different line types correspond to simulations with /t*//? = 0.1 and Ve = 1468 kms~ 
(solid line), k^/k = 0.5 and v^ = 561kms~ (dotted line), k^/k = 0.8 and Ve = 344kms~ 
(dashed line), and k^/k = 0.3 and -Og = 2082 kms~^ (dot-dashed line). 
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Fig. 5. — Relative probability distribution for the physical thin disk scale size r^ for Q0957. 
The solid line represents the probability distribution directly arising from the microlensing 
simulations, while the dashed line shows the result of imposing a prior on the mean microlens 
mass of 0.1 < {M/Mq) < 1.0. 
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Fig. 6. — Acc retion disk size versu s supermassive black hole mass relation (thick solid line) 
and data from lMorgan et al.l (120 lOl ) including the new measurement of r^ for Q0957, scaled to 
2500 A and corrected to 60° inclination. Q0957 is consistent with the mean trend. The dotted 
line shows the scale radius as a function of central black hole mass predicted by theoretical 
thin disk models (for L/ Le = 1/3 and rj = 0.1), while the diagonal crosses indicate the thin 
disk size predicted by the magnification-correct ed luminosity of the diffe rent quasars. We 



f ind th at Q0957 is consistent with the findings of [Morgan et al.l (|2010[ ) and iBlackburne et al. 



( 120111 ) that the optical continuum sizes of quasars measured through microlensing analyses 
are larger than the sizes predicted by thin disk theory for a given black hole mass. The 
microlensing source size for Q0957 is also larger than the luminosity- cons trained thin disk 



size, s imilar to the findings of lPooley et al.l (l2007| ). lMorgan et al.l ( l2010l ) and lMediavilla et al. 
(|201lh . 



